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Abstract 

A method for computing integrals of polynomial functions on compact symmetric 
spaces is given. Those integrals are expressed as sums of functions on symmetric 
groups. 

1 Introduction 

Let G be the unitary group U(?2), orthogonal group 0(n), or (compact) symplectic group 
Sp(2n), equipped with its Haar probability measure, and suppose that G is realized as 
a matrix group. Consider a random matrix X = {xij) picked up from G. Moments of 
matrix elements 

^^iih ■ ■ ■ Xidk^i'.ji ■ ■ ■ XiUl] for G = U(n), 

E[xi^j^ ■ ■ -Xij^jJ for G = 0(n) or Sp(2n) 



have been one of main interesting objects in random matrix theory. Weingarten |Wj 
tried to understand the moments (11.11) for large n. Collins |Coj established a method of 
computations of (II. ip for U{n) with all finite n and called his method the Weingarten 
calculus. The moments are expressed as sums in terms of class functions Wg^(-; n) on the 
symmetric group Sk- The function Wg^(-; n), called the unitary Weingarten function, has 
rich combinatorial structures involving Jucys-Murphy elements |MNt [N] . The Weingarten 
calculus for 0{n) was constructed in |CSnt [CM] , and the corresponding Weingarten func- 
tions, called orthogonal Weingarten functions, are if^-biinvariant functions on S2k, where 
Hk is the hyperoctahedral subgroup of S2k- The orthogonal Weingarten function also 
bears nice combinatorics |MH [Z]. The Weingarten calculus for Sp(2n) was mentioned in 
|CSn|[CSt| . See also a recent preprint [Da], in which Brownian motions on classical groups 
are studied. In the last decade, Weingarten calculus has been widely applied: Harish- 
Chandra-Itzykson-Zuber integrals \Co\ IGGJj . quantum information theory |CN| ICFNj . 
designs [S] , and statistics |CMS ' . 



In his pioneering works for statistical mechanics, Dyson |Dy| introduced three impor- 
tant classes of random unitary matrices as modifications of Gaussian matrix ensembles. 
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These unitary matrix ensembles are known as circular orthogonal/unitary /symplectic en- 
sembles (COE/CUE/CSE). Whereas the CUE is nothing but the unitary group equipped 
with the Haar measure, the COE and CSE are not Lie groups but compact symmetric 
spaces. 

Generally, a classical compact symmetric space is of the form G/K, where G is a 
unitary, orthogonal, or symplectic group, and K is a. closed subgroup of G consisting of 
fixed-points of an involution Q on G. Let S be the image of the map G 3 g ^ Q{g)^^g. 
We thus obtain a matrix realization of G/K: 

G/K^S; G3 g^Q{g)-^g eS, K = {h e G \ Q{h) = h}. 

Moreover, a G-invariant probability measure on S is induced from the Haar measure on 
G via this map. Then the matrix space S equipped with the probability measure is the 
matrix ensemble associated with G/K. In the case of the COE, G/K = U(n)/0(n), 
^{9) = V (the complex-conjugation of g), and S is the set of n x n symmetric unitary 
matrices. Cartan |Ca] classified classical compact symmetric spaces into seven infinite 
series labeled by A I, A 11, A 111, BD I, C 1, C 11, and D 111 (see Figure [1]). Thus we 
obtain seven series of random matrix ensembles associated with these compact symmetric 
spaces. Eigenvalue distributions for them have been studied, see, e.g., [Duj . |CaMj . and 
[F| §3.7.2 and 3.7.3]. In contrast, distributions for those matrix elements were not studied 
so much. Ones can see a small work for only classes A I and A II in [0 §2.3.2]. 

In this article, we will construct Weingarten calculus for random matrices associated 
with each class C = A I, A II,... in Figure [TJ A development in a similar direction can 
be seen in [CStj . but only asymptotic behaviors of two-degree moments were observed 
in that article. As in the case of classical groups U(n), 0(n), Sp(2n), we will express the 
moments of random matrix elements as sums over symmetric groups. To do it, we need 
a distinguishing Weingarten function Wg'' for each C. 

This article is organized as follows. In Section 2, we review the Weingarten calculus 
for U(?7.), 0{n), and Sp(2n), developed in [Co\ \CM\ ICSnt [CSt] . Especially, the case for 
symplectic groups is discussed in detail. It is described in [CStj ID a] , however, unlike their 
descriptions, we introduce the symplectic Weingarten function in terms of twisted zonal 
spherical functions for the twisted Gelfand pair {S2k, Hk,e\H^)- In Section 3, we discuss 
the Weingarten calculus for random matrices from class A I and A II, which are the COE 
and CSE, respectively. The COE case (C = A I) is already given in |M2t IM3] in detail 
but we present it here again for readers' convenience. In Section 4 and 5, we discuss the 
Weingarten calculus for chiral ensembles (A III, BD I, C II) and Bogolioubiov-de Gennes 
(BdG) ensembles (D III, C I), respectively. In the final section, we give a short conclusion 
for new Weingarten functions. 
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Class C 


Symmetric spaces 


Random matrices 


A I 
A II 


U{n)/0{n) 
U(2n)/Sp(2n) 


circular orthogonal ensemble (COE) 
circular symplectic ensemble (CSE) 


A III 
BD I 
C II 


U(n)/U(a) X U(6) {a + b = n) 
0(n)/0(a) X 0(6) {a + b = n) 
Sp(2n)/Sp(2a) X Sp(26) (a + 6 = n) 


chiral ensemble 


D III 
C I 


0{2n)/l]{n) 
Sp(2n)/U(n) 


Bogolioubov-de Gennes (BdG) ensemble 



Figure 1: List of classical compact symmetric spaces. Remark that we consider the 
full groups U{n) and 0(n) rather than special ones SU(r;,) and SO{n), the latter groups 
of which are usually used in Lie theory. 



2 Weingarten calculus for classical groups 
2.1 Weingarten calculus for unitary groups 

In this subsection, we review the Weingarten calculus for unitary groups. See |Col ICSnl 
[CMS] for details. 

2.1.1 Partitions and cycle-types 

A partition A = (Ai, A2, . . . , A;) of a positive integer is a weakly decreasing sequence of 
positive integers with |A| = ^'^^ A, = k. We write £(A) for the length / of A. When A is 
a partition of k, we write A h fc. 

Let Sk be the symmetric group acting on [k] = {1,2, . . . , k}. A permutation a E Sk is 
usually expressed in the two-array notation (^(\) ^^2) o-ffc))' but transpositions are often 
written as {i j) shortly. Denote by id^ the identity permutation in Sk- If a permutation a 
in Sk is decomposed into disjoint cycles with lengths /ii > /i2 > • ■ ■ > /W/, then the partition 
/i = (/ii,/i2, . . . , /i;) of /c is called the cycle-type of a. For example, the cycle-type of the 
permutation ( 3 f 2 I e 5 ) in S'g is (3, 2, 1). 

2.1.2 Unitary Weingarten functions 

Let L{Sk) be the algebra of complex-valued functions on Sk with convolution 
(/i * /2)(^) = Yl M^)Mr-'^) (/i, /2 e L{Sk), a G Sk). 

The identity element in the algebra L{Sk) is the Dirac function did,.- 

Let Z{L{Sk)) be the center of L(^fc): Z{L{Sk)) = {h e L{Sk) \ h * f = f * h {f E 
L{Sk))}. For a complex number z, we define the element T^(-; z) in Z{L{Sk)) by 



TU(a;z) = /('^) {aeSk), 
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where /i is the cycle-type of a. Note that the upper index U stands for the unitary group. 
The class function T^(-; z) can be expanded in terms of irreducible characters of Sk as 
follows: 

Ahfc 

where := x^iy^k) is the dimension of the irreducible representation associated with A 
and C\{z) is the polynomial in z given by 

Cx{z)= n + 

Here G A stands for 1 < z < £(A), 1 < j < Aj. In other words, the (i, j) are 

coordinates of the Young diagram of A. 

The unitary Weingarten function Wg^(-; z) on Sk is defined by 

Ahfc 
Ca(z)7^0 

summed over all partitions A of A; with C\{z) 7^ 0. It is the pseudo- inverse element of 
;z), i.e., the unique element in Z{L{Sk)) satisfying 

TU(-;;^)*WgU(-;z)*TU(-;z)=TU(.;z). 

In particular, unless z G {0, ±1, ±2, . . . , ±(A; — 1)}, functions T^(-; ,2) and Wg^(-; z) are 
inverse of each other and satisfy T^(-; z) * Wg^(-; z) = 5idj.. 



2.1.3 Integrals on unitary groups 

The unitary group U(n) is \i{n) = {U G GL(r2,C) | UU* = In}, which has the Haar 
probability measure. Here U* := is the adjoint matrix of U and = (5ij)i<ij<n 
is the n X n identity matrix. We will also write U for a random element of U(n), with 
distribution given by Haar measure. 

Theorem 2.1 (Weingarten calculus for unitary groups |Cot ICSn] ). Let U = {uij)i<i,j<n 

be an n X n Haar unitary matrix. For four sequences i = {ii, . . . ,ik), j = (ji, . . . ,jk), 
i' = {i'l, . . . ,i'k), j' = Ui, ■ ■ ■ ,jk) of positive integers in [n], we have 

Here 6a{i,i') is defined by 

k 

<5.(i,0 = n'^M.).^^- (2-2) 

s=l 
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If k ^ I, then E[uij^j^ ■ ■ ■ Ui^^j^Ui'^f^ ■ ■ -Wij^-'] vanishes for any indices ii, . . . ,ik, ji, ■ ■ ■ ,jk, 

^1? • • • 5 ^i; Jl5 • • • ) J/ ■ 

Example 2.1. Wg^(idi;n) = i; 

WgU(id2; n) = -; WgU((l 2); n) ~^ 



(n + l)(n-l)' " ' n(n+l)(n-l)' 

We can see more examples in |Co] . 



Remark 2.1. In the sum of (12.11) . we can drop the restriction for A. Even if we replace 
the definition by 

which is a rational function in ^ G C with finitely many poles, then Theorem 12.11 remains 
true after cancellations of poles. This enables us to ignore the restriction C\{z) 7^ in 
latter sections. See |CSnt Proposition 2.5] for details. 



2.2 Weingarten calculus for orthogonal groups 

In this subsection, we review the Weingarten calculus for orthogonal groups, developed 
in |CSnt ICMt ICMSj . The theory of zonal spherical functions for finite Gelfand pairs are 
seen in |Mac[ Chapter VII]. 



2.2.1 Hyperoctahedral groups and coset-types 

Let Hk be the hyperoctahedral subgroup of 5*2^ with order 2^k\, generated by adjacent 
transpositions {2i — 1 2i) {1 < i < k) and double transpositions {2i — 1 2j — l)(2i 2j) 
{1 < i < j < k). Let M2k be the subset of permutations cr G 5*2^ satisfying 

a{2i - 1) < a{2i) {1 < i < k) and (t(1) < a(3) < ■ ■ ■ < (t(2A; - 1). 

The elements a in M2k form the complete set of the representatives of S2k/Hk and are 
sometimes identified with perfect matchings of the forms 

{{a(l), a(2)}, {^(3), a(4)}, . . . , {a{2k - 1), a{2k)}} 

on [2k]. 

For a permutation a G S2k, we define the graph r(cr) as follows. The vertex set is 
{1, 2, . . . , 2k} and the edge set consists of red edges {2r — 1, 2r} and blue edges {a{2r — 
l),(T(2r)}, where r runs over l,2,...,k. Then all connected components of the graph 
r(o") have even vertices with numbers 2/ii > 2/i2 > ■ ■ ■ > 2fii. We call the partition fi : = 
(/ii,/i2, ■ ■ ■ , f^i) t~ k the coset-type of a. For example, the coset-type of the permutation 
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(312465) 'S'e is (2, 1). The coset-type distinguishes with double cosets of Hk in S2k- 
Specifically, for two permutations a, r in S2k, their coset-types coincides if and only if 
HuaHu = HkTHu VII, (2.1)]). 

For a partition /i = (/ii, /i2, . . . , /i;) of fc, we define the permutation in 5*2^ as follows: 
For each r = 1, 2, . . . , /, 

(2 Eg + 2) = 2 EJli + 2/i., (2.3) 

(2 EI=i Ati + P) =2 EI=i yUi + P - 1 for p = 3, 4, . . . , 2/i^. 

For example, cr(3 1) = (i6234578)- -'■^i^ ^^^y ^^^^ permutation cr^ belongs to 

and has the coset-type /x and signature +1. In particular, cT(ifc) = id2fc, the identity 
permutation in S2k- We often regard as a typical permutation of coset-type /x. 

2.2.2 Orthogonal Weingarten functions and Gelfand pairs 

Let L{S2k,Hk) be the subspace of all iffc-biinvariant functions in L{S2k)'- 

L{S2k,Hk) = {fe L{S2k) I fiCcr) = f{aO = f{a) (a G ^2^, C e H^)}. 

It is well known that {S2k, H^) is a Gelfand pair, i.e., L{S2k, H^) is a commutative algebra 
under the convolution f jMacl VII, (1.1) and (2.2)]). 

For two functions /i, /2 in L{S2k), we define a new product /i t^t /2 by 

(/i ^ /2)(a) = J2 hi^)h{r-'a) (a G ^2^). 
Example 2.2. For /i,/2 G ^(^4), 



(/l 




/2)( 


1 
1 


2 3 4 \ 
2 3 4 J 


= /ll 


^ 1 


2 3 4 
2 3 4 


)/2( 


' 1 2 3 4 ^ 
. 1 2 3 4 J 




' 1 


iil)/2i 


2 3 4 ^ 
^ 1 3 2 4 J 


1 + 


^ 1 


lit)/2( 


' 1 


2 3 4 ^ 

3 4 2 J 


(/l 




/2)( 


1 
1 


2 3 4 N 

3 2 4 J 


= /ll 


^ 1 


2 3 4 
2 3 4 


)/2( 


' 1 2 3 4 ^ 
. 1 3 2 4 J 




^ 1 


3 2 4 ) /2 1 


^ 1 2 3 4 ^ 
2 3 4 J 




' 1 


4 2 3 ) /2 I 


' 1 


2 3 4 N 
43 2 J 


(/l 




/2)( 


1 
1 


2 3 4 N 
4 2 3 J 


= /ll 


' 1 


2 3 4 
2 3 4 


)/2( 


' 1 2 3 4 ^ 
. 1 4 2 3 J 




' 1 


iil)/2i 


2 3 4 ^ 
^ 1 4 3 2 J 


1 + 


' 1 


lii)/2( 


' 1 


2 3 4 ^ 
2 3 4 J 



For /i,/2 G L{S2k,H,), we have /i ^ /2 = (2'=A;!)-7i * /2 G L{S2k,Hk). In fact, 
since M2fc is the complete set of representatives of cosets aHk in S2k and since /i, f2 are 
iffc-biinvariant, we have 



Hence L{S2k, Hk) is a commutative algebra under the product -k with the identity element 




1 if a G Hk, 

ifaG52A^fc- 
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For each partition A = (Ai, A2, . . . ) l~ k, we define the zonal spherical function of 
the Gelfand pair {S2k, Hk) by 

with 2A = (2Ai, 2A2, . . . ). The u^, X\- k, form a hnear basis of L{S2k, Hk) and satisfy the 
orthogonahty relation 

-'*-''-^J^^y (2-4) 

(IHaa VII, (1.4)]). 

For a complex number z, we define the element T'-'(cr; z) in L{S2k, Hk) by 

T°(a;^) = /('^) iaeS^k), 

where n is the coset-type of a. We emphasize that T° is different from T^. The function 
T*-*(-; z) is expanded in terms of u'^ as follows ( |CMt (4.5)]): 

\\~k 



where Dx{z) is the polynomial in z given by 



Dx{z)= n {z + 2j-i-l). 



The orthogonal Weingarten function Wg'^(-; z) on S'2fc is defined by 



(2fc)! ^ D,[z, 



which is the pseudo-inverse element of T'^{-]z), i.e., the unique element in L{Sk,Hk) 
satisfying 

T0(-; z) * Wg°(-; z) *T0(-; z) = T0(-; z). (2.5) 

In particular, if D\[z) ^ for all partitions A of A;, functions T°(-; z) and Wg'~'(-; z) are 
their inverse of each other and satisfy T*-'(-;z) *Wg*~'(-;2;) = 1^. Those relations follow 
from (12. 4p and expansions of T°(-; z) and Wg°(-; z) in terms of u^. 

2.2.3 Integrals on orthogonal groups 

The (real) orthogonal group is 0{n) = {R G GL(r?,,]R) | RR^ = /„} and has the Haar 
probability measure. 
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Theorem 2.2 (Weingarten calculus for orthogonal groups |CSnllCM] ). Let R = {rij)i<ij<n 

be an n X n Haar orthogonal matrix. For two sequences i = {ii, . . . ,i2k) and j = 
(ii; • • • ) J2fc) of positive integers in [n], we have 

nr^^n■■■r^,,JJ= J2 A<.(*)Ar(j)WgO(a-V; n). 

a,T€M2k 

Here Ao-(*) is defined by 

k 

A.(i) = nWi)=M..)- (2.6) 

3 = 1 

Furthermore, IE[rjj^ ■ ■ ■ rj^^+ijafe+il = M any ii, . . . , i2k+i,ji, ■ ■ • , J2fc+i- 
Example 2.3. Wg°(id2;n) = 



Wg^(id4;n) = . , -, Wg"(a(2);n) 



n{n + 2){n-iy ^ ' ^ " ' n(n + 2)(n-l) 

We can see more examples in |CSn[ ICMj . 



2.3 Weingarten calculus for symplectic groups 

In this subsection, we give the Weingarten calculus for symplectic groups. It was given 
in |CStt ID a] , however, unlike their descriptions, we employ the theory of twisted Gelfand 
pairs ([Macl VII, Examples 1-10, 1-11, 2-6, and 2-7]). 

2.3.1 Twisted Gelfand pairs 

Let e be the signature function on S^k and consider the linear space 

U{S2k. Hk) = {fe L{S2k) I /(C^) = /«) = e(C)/(^) G S2k, C e Hk)}. 

The space L''{S2k, Hk) is closed under the convolution *, and it becomes a C-algebra. 
Furthermore, it is known that L'^{S2k, Hk) is a commutative algebra, which means that the 
triple {S2k, Hk, el^j,) is, by definition, a twisted Gelfand pair ( |Mact VII, Example 2-6]). It 
is immediate to see that, if /i, /2 G L^(S'2fc, Hk), then /i*/2 = {2^k\)^^ fi* f2 G L''{S2k, Hk). 
Thus, L^{S2k, Hk) is a commutative algebra under the product ★ with the identity element 



m^) 



e(o-) if a G Hk, 

if a G 52fc \ Hk. 



For each partition A = (Ai, A2, . . . ) l~ fc, we define the twisted spherical function tt^ of 
the twisted Gelfand pair by 

= i2'^k\)-'x''" * 11 
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with A U A = (Ai, Ai, A2, A2, . . . ). The map 

L{S2k, Hu) -> U{S2k, Hk) -.f^r, rW) = e{a)f{a) {a G S2k) 

defines a C-algebra isomorphism. The twisted spherical function vr^ is the image of u^' 
with vr^((T^) = (cr^)- Here A' = (A'^, Aj, . . . ) is the conjugate partition of A, which is 
characterized by G A <^ G A'. Hence, the {vr^ \ X \- k} form a hnear basis of 
L^{S2k, Hk) and satisfy the orthogonahty relation 



(2fc)! 1 ,A 



Tr^^TT^" = 6x,,--—-——TT 



2.3.2 Symplectic Weingarten functions 

Let z be a complex number and consider the function T^p(-; z) in U{S2k, Hk) defined by 

TSp(a; z) = (-l)^e(a)(-2z)^(^) {a G ^2^), 
where is the coset-type of a. Let 

{«,i)GA 

The expansion of T^^(a; z) in terms of linear basis vr^ is given as follows. 
Lemma 2.3. 

^''^■'^) = ^E/'''^AW-^ (2.7) 

Proof. We first suppose that z = n, a positive integer. We employ symmetric function 
theory. Recall power-sum symmetric functions and twisted zonal functions Z'^. We only 
use the following properties for them here (see |Mact VII, Example 2-7]). For partitions 
A, /i of fc. 



Xhk 



Here (1") = (1, 1, . . . , 1) with n times. Hence, if fi is the coset-type of cr G S2k, then 



T^-{a-n) = e{a){-l f-'(^h'^^'>p,{r) = ^(a)^ /^^V(a,)Dl(n), 

AhA: 



which implies the desired formula for 2; = n. 

Since the both sides on (12. 7p are polynomials in z, the equalities at all positive integers 
z = n implies the ones at all complex numbers z. □ 
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The symplectic Weingarten function with parameter z is the function in U{S2kiHk) 
defined by 

ofc^l , fAUA 

Note that rSp(a;z) = {-l fe{a)T^{a--2z) and Wg^P(a;z) = (-l)^e(a)Wg°(a; -2^). 
Equation (12.51) is equivalent to 

TSP(-;;z)^Wg^P(-;^)*TSP(-;^) = T^''(-;-2)- (2.9) 

2.3.3 Integrals on symplectic groups 

Consider the vector space C^" of column vectors with standard basis (ei, 62, . . . , e2n). 
Define the skew-symmetric bilinear form (-, ■) on by 

{v,w) = v^Jw (2.10) 

with 

Let (e^, 62 , • • • , e2„) be the dual basis of (ei, . . . , e2n) with respect to (■, ■): 

(er,e,)=5,, {i.je[2n]). 

More specifically, (e^, . . . , e^, e^.^^, . . . , e^„) = (-e„+i, . . . , -Cgn, ei, . . . , e„). For conve- 
nience, we use the following notation: 

{1 if 1 < 2 < n and j = i + n, 
-1 if 1 < j <nandi = j + n, 
otherwise. 

Note that 

2n 2n 

e, = J^et = Y,{i,v)el, = Je, = ^(p, i)ep (^ = 1, 2, . . . , 2n) (2.11) 

and J = ((«,j))i<ij<2„. 

For a 2n X 2n matrix X, we define the dual matrix of X by 

X° := JX^JT. 

Then = {X^v,w) for all G C^^^. 
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We realize the (unitary) symplectic group Sp(2n) by 

Sp(2n) = {S= {s,,)i<,,,<2n G U(2n) I 55° = hn}. 

It is equipped with the Haar probabihty measure. The following theorem was first given 
in [CSt] without the explicit expression ( 12. 8 P for Wg^^. 

Theorem 2.4. Let S = (■Sij)i<ij<2n be a Haar symplectic matrix. For two sequences 
i = (ii, . . . , i2k) and j = {ji, ■ ■ ■ ,j2k) of positive integers in [2n], we have 

o;T€M2k 

Here the symbol A'^{i) G {0, 1, —1} is defined by 

k 

AU*):=n<M2r-i),M20)- (2-12) 

r=l 

Furthermore, E[si^j^ ■ ■ ■ Sj^^+ja^+J = for any ii, . . . , ^2fc+l, Ji, • • • , J2fc+i- 

We postpone the proof of this theorem in the next subsubsection. 
Example 2.4. Wg^P((T; n) = ^ for a G ^2; 

2n — 1 

Wg^P(a; n) =e(a)— for a E H^; 

4n[n — l)[2n + 1) 

WgSP(cr; n) =e{a)— \ for a e S,\ H^. 

4n(n — lj(2n + I) 

2.3 A Proof of Theorem Q 

This proof was given by Collins and Stolz |CSt] (see also |CM] ). however, they did not 
give any explicit expression for the symplectic Weingarten function. We here reconstruct 
their proof and observe how Wg^'' arises. 

We first introduce the following useful notation. For a permutation a G and two 
matrices X = {xij)i<ij<n and Y = {yij)i<i,j<n, put 



k k 

(2r-l)''J(T-l(2r) ' 



T^Ct{X.,Y) — jj^ ^Pa{2r-l).P<T(2r)?/p2r-l,P2r ~ ^l?2r- 1 ,92,- 1/?^- 1 

Piv,P2fcG[n] r=l gi,---,g2fc6M ^'=1 



Note that %-i (X, Y) = %{Y, X) and (X, Y) = HS (X, Y), where is defined 



in (12. 3p . Since 

'^(m) ^) ~ ^ ] '^Pl,P2ml/pi,P2'^P2,P3l/p3,P4 ' ' ' '^P2m-2 ,P2m- 1 ' ' ' l/p2m-l,P2m5 



Plv,P2m 
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we have 

T {X Y) I '^^K"^^)'"] if X is a symmetric matrix, 

1 — Tr[(Xy)"*] if X is a skew-symmetric matrix. 

As a particular case, we can see 

r.(J, J) =TSP((7;n) {aeS2k). (2.14) 

Indeed, the skew-symmetry of J imphes that the function a i— )• Ta{J, J) belongs to 
L^{S2ki Hk), and hence it is enough to check Ta^{J,J) = T^^{afj,]n) for all partitions /i. 
However, it follows from (EH that {J, J) = T.,,^, ( J, J) = IlS [" Tr(-/2„)'^'] = 
(_l)fc-^(M)(2n)^(^) = TSp(a^;n). 

We next recall the invariant theory for symplectic groups. Consider the tensor product 
^£2n-^(i)2k g^^^ define a bilinear form on (C^"^)®^'' by 

2k 2k \ 2k 

'^Vj,(^Wj \ := {vi,...,V2k,Wi,...,W2k G C^"), 

,i=l j=l I j=l 

where the skew-symmetric bilinear form {v,w) on the right hand side is defined in (12.101) . 
Note that this bilinear form on (C^")'^^'^ is symmetric. 
Put 

9k = Yl < ® ® ■ ■ ■ ® ® e,, e {e-r^\ 

pi,...,Pk&[2n] 

The symmetric group S2k acts on (C2n^(»2fc 
while the symplectic group Sp(2n) acts by 

S{Vi®---® V2k) = Svi (g) ■ ■ ■ Sv2k- 

Then the First Fundamental Theorem for the symplectic group states that {p2k{(^)0k \ o" G 
M2k} spans the vector subspace of (C^")®^'^ consisting of invariant elements under the 
action of Sp(2n) (see, e.g., [GWl Theorem 5.3.4]). 

Let us go back to the proof of Theorem 12. 4[ Let G be the symmetric matrix 

G = {{p2k{cr)0k, P2k{r)9k))a,TeM2k (2.15) 

and W = (w((T, r))o-_T-gAf2fc the pseudo-inverse matrix of G, i.e., W is the unique symmetric 
matrix satisfying GWG = G. 



Weingarten calculus for compact symmetric spaces 



13 



Let ^ = f Sij)i<i,j<2n be a Haar symplectic matrix. Since each matrix element Sij is 
expressed as Sij = (e^, Scj), we have 

Ehiii ■ ■ ■ Si2kj2k] =E (e,^, ® ■ ■ ■ ® el^,S{ej, ® ■ ■ ■ O Cj^J) 

= (< ® ■ ■ • ® <,,E[5(e,, ® ■ ■ ■ ® e,,J]). 

As mentioned above, for any v G V^^'', the Sp(2n)-invariant vector E[S'f ] can be expanded 
in terms of p2k{o')0k- From a discussion parallel to the proof of |CMt Theorem 2.1] (see 
also |CSt] ). the expansion is given by 

¥.[Sv] = ( Yl w(a,r)(^;,p2fc(r)^fc) J P2fc(a)^fe. 

Applying this to the previous equation, we obtain 

cj,TeM2k 

Here the values of bilinear forms are computed as follows: 

{el ® • ■ ■ ® el^,p2k{<y)ek) = {P2k{<y-^){el ® ■ ■ ■ ® e.^^J, 9^) 

k / 2n ^ 
r=l \p=l 

and, similarly, (e^^ ® ■■• ^ ej^^,, p2k{T)9k) = A'^(j). In conclusion. Theorem 12.41 follows 
from the next lemma. 

Lemma 2.5. For a,T e M2k, w(a-, r) = Wg^P(a"V; n). 
Proof. If we put 

T(cr) = {9k, p2k{o)9k) (o- G 52^), 
then G = (T(cr^^r))o-,rGAf2fc- Using fl2.1ip we have 

^k= Y (Pl'P2)(P3,P4) ■ ■ ■ (P2fc-l,P2fc)ep^ ® (g) ■ ■ ■ (g) e. 



V 

'P2k 



Pl,---,P2k 



Y ^^1' ^2) (93, 94) ■ ■ ■ (g2fc-i, g2fe)egi (8) o ■ ■ ■ ® 



'?2fe' 



and hence 

T{(y) = Y (^'^{i)'^'^{2)) • ■ ■ (9a(2fc-i),g<7(2fc))(gi,g2) • ■ ■ {q2k~i,q2k) = %{J, J), 

'?lv,'?2fe 

which implies T((t) = T^p(o"; n) by (12.141) . The matrix (Wg^^(cr^^r; n))^^reM2k is therefore 
the pseudo-inverse matrix of G = (T^P((T~^r; n))a,T&M2k by (12. 9p . □ 
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3 Circular ensembles 

From now we consider random matrix ensembles S associated with classical symmetric 
spaces G/K. As we mentioned in Introduction, such ensembles are realized in the following 
way. 

G/K ~ 5; G3 ^{g^^g G S. 

Here Q is an involution on G and K is the fixed-point set of Q. If X is a Haar random 
matrix picked up from G, then the matrix V := Q{X)~^X is a random matrix associated 
with G/K. We consider the seven series of random matrices associated with compact 
symmetric spaces given in Figure [H 

In this section, we deal with the most important classes: circular orthogonal ensembles 
(COE) and circular symplectic ensembles (CSE). 

3.1 Class A I 

The setting for A L G = U{n), K = 0{n), Q{g) = 'g. S consists of n x n symmetric 
unitary matrices. 

When U is an n x n Haar unitary matrix, a random matrix corresponding to U(n) /0{n) 
is defined by 

V = V^^ = VL{U)-^U = U^U 

and is said to be a COE matrix. The Weingarten calculus for a COE matrix is constructed 
in |M2j . For completeness of this paper, we review main results in |M2j . Applying the 
Weingarten calculus for U(n), we can obtain the following theorem. 

Theorem 3.1 (Theorem 1.1 and Proposition 3.1 in |M2j ). Let V = V^^ = ivij)i<i,j<n be 
an n X n COE matrix. For two sequences i = {ii, . . . , i2k) and j = (ji, . . . ,j2k), we have 

with the convolution Wg^^{-;n) := T^{-;n) * Wg^(-;n) in L{S2k)- If k ^ I then 



^['^ili2 '^isM ■ ■ ■ ^i2k-li2k^jlj2^33H ' ' ' ^321-13211 

always vanishes. Moreover, the function Wg^^{-;n) belongs to L{S2k,Hk) and coincides 
with the orthogonal Weingarten function with parameter n + 1, i.e., 

Wg^'(a;n) = Wg°(-;n + l). 
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As a corollary of this theorem, the following moments for a single entry are computed 
(see [M2l Theorems 4.1 and 4.2]): 



'{n + l){n + 3)---{n + 2k-r 



^f'""'^''' ] ^n{n + l){n + 2)---{n + k-2){n + 2k-l) ^'^^^^ 

The first equation can be obtained easily from the previous theorem, but the derivation 
of the second one is somewhat complicated. 

Example 3.1. From Example 1X51 we have Wg^^{a] n) = Wg^{a; n+1) = ^ for a G 5*2; 



n+2 



for a E H2 



Wg^\cr; n) = WgO(a; n + 1) = { -("+i)("+3) 

. n(n+l)(n+3) for (7 G ^4 \ 



3.2 Class A II 

The setting for A II: G = U(2ra), K = Sp(2n), n{g) = {g^y^. S consists of 2n x 2n 
unitary matrices g satisfying g^ = g. 

When f/ is a 2r;,x2?T, Haar unitary matrix, a random matrix corresponding to U(2n)/Sp(2r;,) 
is defined by 

V = V^'' = U^U 

and is said to be a CSE matrix. 

We would like to compute mixed moments for Vij and vij. In order to simplify the 
notation, we deal with 

Vij := (ci, Vcj) (1 <i,j < 2n) 
instead of Vij = (e^, Vcj). More specifically. 



Vij 



-Vi+n,j ifl<i<n, 



V 



t—n 



J if n + 1 < i < 2n. 



Theorem 3.2. Let V be a 2n x 2n CSE matrix. For two sequences i = {ii, . . . ,i2k) (md 
3 = (ji, • • • , J2fe) in [2n]''^'', we have 

E[^n,i2^i3,i4 ■ ■ ■ Vi,,_„i,,Vj,j,Vj^j, ■ ■ ■ Vj^,^„jJ = J2 j)Wg^"((T;n) 



wt/i Wg^"(-;n) :=TSp(.;n) * WgU(-;2n). If k I then 
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Proof. Each element of a 2n x 2n CSE matrix V = U^U is given as 

Vi^i' = {ei,U^Uei>) = {Uei,Uee) = ^ UpiUgi'{p,q) {i,i'e[2n]). 

p,qe[2n] 

Hence 



^Nl,i2'^j3,i4 ■ ■ ■ ^«2fe-l,«2fc'^il J2'^j3j4 ' ' ' '^J2fc-1 ,i2fc] 

p=(pi,...,P2fc) q={giv,g2fc) ''=1 

(When k ^ I, WyUp^i^ ■ ■ ■ Up^,iJ^^J^'^^'^u^^ = 0.) Applying the Weingarten calculus for a 
Haar unitary matrix U, we have 

E[5n,i2^i3,M ■ ■ ■ '^i2.-i,i2fe^ii,i2^'j3,j4 " " " ^^^2.-1^2.] = ^^^^'^^ Y Wg^ (t" V; 2r2)f (t) , 
where 

fe 

P=(pi,-,P2fc)G[2n]x2fc q=(gi,...,q2ft)g[2n]x2fc r=l 

=r.(J,J) = rSp(r;n). 

The last equality above follows from (12 .Mp . □ 

Proposition 3.3. The function Wg^^^{-;n) = T^P(-;n) * Wg^(-;2r;,) coincides with the 
symplectic Weingarten function with parameter n — ^. Specifically, for each a G 5*2^, we 
have 

Wg'^"((T;n) = Wg^P((T— i~ 
Proof It follows from (EZD and (El]) that 



Ahfc ^ ' /xh2fc 



Since tx^ * = Sxux,f,^^T<'^ and since 
/^U^) _ n(„-)eA(2^-2^+J + l) 



2^* 



Caua(2^) n(,,)eA(2^ + J - (2^ - 1))(2^ + J - 2^ n(„-)eA(2^ + J - 2^ " 

we have T^^{-, n) * WgU(-; 2n) = Eak. I^fS^^^' = Wg''( ' ^ - ° 

Example 3.2. From Example [23] we have Wg'^"((j;n) = Wg^P(cT;n - |) = for 
cr G So; 



Wg^"(a;n) = Wg^P(a;n-i; 



<^) 2n(2n-l)(2n-3) for a G ^4 \ i^s- 
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4 Chiral ensembles 

In this section, we deal with random matrix ensembles associated with classes A III, BD I, 
C II. They are known as chiral ensembles. 

4.1 Class A III 

The setting for A III: G = U{n), K = U(a) x U(6), n{g) = I'^f^gl'^f,, where n = a + b with 
a > 6 > 1 and 

'la O 

O -h, 

Let U he annxn Haar unitary matrix. A random matrix corresponding to U(n) /U(a) x 
U(6) = SU(n)/S(U(a) x U(6)) is defined by V = 1^^™ = r^t,U*r^^U. For the sake of ease, 
we consider a Hermitian and unitary random matrix 

w = W^^^^ = U*I'^^U, 

instead of = I'^^W. 



T' 



We define the function Tj|"i in Z{L{Sk)) by 



where Tr<,(A) = Uf^l Tr(A^O if 

/i is the cycle-type of a. More specifically, we have 
T„Vn(a) = (a + 6f(^)(a-6f(^). 
Here (resp. ^°(/i)) is the number of parts fij with even lengths (resp. odd lengths). 

Theorem 4.1. Let W = W^^^^ be the random matrix defined as above. For two sequences 
i = andj = {ji,...,jk), we have 

E[wi,j,Wi^j^ ■ --Wi^j^] = j)Wg^"^(a; a, b). 

Here the function Wg'^"^(-; a, 6) in Z{L{Sk)) is defined by 

Wg^"i(-;a,6)=T,V"*WgU(-;n). 

Proof. The random matrix W is Hermitian and the distribution of W is invariant under 
the unitary transform W i— )■ gWg*, where g is a fixed unitary matrix. Hence we can apply 
Theorem 3.1 in |CMSj and obtain 

nWhnWi2h---Wi,j,]= 5.(*,j)WgU(r-V;n)E[Tr.(iy)]. 

Here, wee see that E[Tr^(Vr)] = E[Tr^(/;;,)] = Tr^(/;j = Tj^"^(r), and we obtain the 
desired identity. □ 
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Example 4.1. n = a + b. Wg^™(idi; a, 6) = 



.Aiii/,i.„n {a-b + l){a-b-l) aiii//. o^. „ n ^aft 



(ra + l)(n — 1) n{n — lj[n + 1) 

Recall power symmetric functions and Schur functions sx. They have the relation 

if /i h /c is the cycle-type of a ( [Macj I, (7.8)]). Furthermore, it is easy to see that 

(-1)^) = p,{i, 1, . . . , 1,-1, -1, . . . , -1) = T,V"(o-). 

b 



Hence we obtain the expansion of T^^^^ in terms of irreducible characters x^'- 



On the other hand, from the well-known identity sa(1") = ^—jr^, the unitary Weingarten 



fc! ' 

function is expressed as 

fX\2 

,X 



Consequently, using the relation x'^ * = jr^XfiX^i we obtain the following expansion of 
Wg^™(-; a,b) in terms of x^'- 

asaIi-.I-i)"),, 



Wg""(.; a. b) = T*" . Wg"(.; „) = i f J > X 

Ahfc ^^^^ ' 



4.2 Class BD I 



The setting for BD I: G = 0(n), K = 0(a) x 0(6), n{g) = r^b9Kb^ where n = a + b with 
a > 6 > 1. 

The discussion is parallel to that in the previous subsection. We deal with the sym- 
metric and orthogonal random matrix 



where i? is an n x n Haar orthogonal matrix. This random matrix is associated with the 
symmetric space 0(n)/0(a) x 0(6) = SO(n)/S(0(a) x 0(6)). 
We define the function T^°^ in L{S2k, Hk) by 



T^°V) = r.(/:„/„). 
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More specifically, if a G S2k has the coset-type fi = (/ii, fi2, ■ ■ ■ , fJ'i), then 
TT\^) = n Tr(4)- = (a + br^\a - bf^^l 

i=l 

Theorem 4.2. Let W = W^^^ be the random matrix defined as above. For a sequence 
i = {ii, . . . , i2k), we have 

E[wi,i^Wi^i^ ■ ■ ■ Wi^^_^i^^] = A^{i)Wg^^\cT;a,b). 

Here the function Wg^^^(-; a, b) in L{S2k, H^) is defined by 

Wg^°i(-;a,6)=Tr^^WgO(-;n). 
Proof. The proof is the same with that of Theorem 14.11 Apply Theorem 3.3 in |CMSj . □ 
Example 4.2. n = a + b. Wg^°^(id2; a, b) = 

Wg [[i23,),^,b)- ^^^2)(n-l) ' 
Wg---(( 1 i i I ) ; a, 5) ^ WgB^^(( 1 1 i I) ; a, 6) ^ ^^^ /^^^^ _ . 

Using zonal polynomials Z\ (see [Mad VII-2]), we can obtain the expansion of Wg^°^(-; ^ 
in terms of zonal spherical functions u^: 



Ahfc 



4.3 Class C II 



The setting for C II: G = Sp(2n), K = Sp(2a) x Sp(26), n{g) = I'^^glab^ where n = a + b 
with a > 6 > 1 and 

J" _ i ^ab ^ 



ab 



O I'ab 



The discussion is parallel to the previous subsection again. We deal with the random 
matrix 

W = W^^^ = S^I'^^S. 

where S* is a 2n x 2n Haar symplectic matrix. It is associated with the symmetric space 
Sp(2n)/Sp(2a) X Sp(26). 
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We define the function T^^^ in U{S2k, Hk) by 
If (J G lias tlie coset-type fi, tlien 

i=l i=l 

1=1 

Here tfie second equality above follows by f l2.13p . 

As in the case of class A II, we consider wtj = (cj, Wcj) instead of matrix elements 

Wij. 

Theorem 4.3. Let W = W'"^^ be the random matrix defined as above. For a sequence 
i = {ii, . . . , i2k), we have 

creAhk 

Here the function Wg*"^^(-; a, b) in L^{S2k, Hk) is defined by 

WgC"(-;a,6)=TC"^WgSP(-;n). 

Proof. Since 

2n 

p,q=l 

we have 

k 

Pl,--;P2k r=l 

with {p,q)" := {ep,I"i,eq). The Weingarten calculus for symplectic groups gives 
where 

k k 

T{t) := JJ(p2r-l,P2r)" = YliPr{2r~l) , Pr{2r)) {P2r~l, P2r)" 

p=(pi,...,P2fc) ^'=1 Plv,P2fe '-=1 

□ 
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Example 4.3. n = a + b. Wg^"(id2; a, 6) = 

(-l)Wg^"(( } 1 i 1 ) ; a, 6) = Wg^"(( II i i ) ; a, 6) = ^(^ _ ^j^^^n + 1) " 

Using twisted zonal polynomials Z'^ (see |Mac[ Vll, Example 2-7]), we obtain the 
expansion of Wg^"(-; a, b) in twisted spherical functions tt^: 



5 BdG ensembles 

In this section, we deal with matrix ensembles of types D III and C I. They are called 
Bogolioubov-de Gennes (BdG) ensembles. 

5.1 Class D III 

The setting for D III: G = 0(2n), K = 0(2n) n Sp(2n) ~ U(n), n{g) = (g^^K 
Let i? be a 2n X 2ri Haar orthogonal matrix. We consider the random matrix 

associated to the symmetric space 0(2n)/U(n). 
We define the function T°"^ on 5*2^ by 

It satisfies T^^^^CO = eiOT^^^K^) for any a G S2k and CC e Hk. Note that 

rr"K)=n^r"K..,)=n^'(^'")- ^^^^ 



i=l j=l 



otherwise. 



Here a partition yU is said to be even \i [i = 2v for some partition v. In particular, 
TP™ (a) = (a G 52fc) if is odd. 

Theorem 5.1. Let V = be the random matrix defined as above. For a sequence 

i = (*i, "^2, • • • , i2k) with an even number k, we have 

^KM<M---K,-„i2u]= A,(i)Wg°"i(or;n). 
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Here Wg°™(-;n) is the function on S2k defined by 

Wg°™(-; n) = * Wg«(-; 2n). 

If k is odd, then, for any sequence i = (zi, Z2, • • • , iik), 

Proof. The proof is similar to that of Theorem 14.31 The Weingarten calculus for a Haar 
orthogonal matrix R gives 

where 

k 

fir) = ^r{p) n W-i,P2.) = J) = T^''\t) 

p=(pi,...,P2fc) ''=1 

for any r G S2k- □ 
Note that 

Wg°"^(CaC'; n) = e(C)Wg°™(a; n) (a g S^,, C, C e 

and hence Wg°™(C;Ti) = for C G i^fc- Furthermore, Wg°™((T;ra) = (ct G ^2fc) if k is 
odd. 

Example 5.1. Since = {3 4) {HH) and {3 A) e H2, we have 

-1 



(-l)Wg^"^((liil);n)=Wg^"^((lli^; 



n] 



2n-l 
5.2 Class C I 

The setting for C I: G = Sp(2n), K = U(n), Q{g) = I'^^gl'^^. Here the fix-point set of Q 
in G is 

{(o f) I ^eU(n)}^UH. 
Let S" be a 2n X 2n Haar symplectic matrix. We consider the random matrix 

instead of V*^^ = Q{S)^^S = I'^^^S^I'^^S, which is associated to the symmetric space 
Sp(2n)/U(n). 
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We define the function ^ on S2k by 

Since J/^„ = ~ ( /" o" ) symmetric and J is skew-symmetric, the function T^^ satisfies 
'(CO = 7;Ci(a)e(C') for (7 e S2k and C, C e Hk. For a = a^, we have 

i=l 2=1 i=l 

/i IS even, 
otherwise. 

In particular, T^\a) = r„D"i((7-i) = ±T^^^\a) for all cr e Sik- 

Theorem 5.2. Le^ = W^^ be the random matrix defined as above. For a sequence 
i — iii,i2, ■ ■ ■ , i2k) with an even number k, we have 

Here Wg*^^(-;n) is the function on defined by 

WgCi(-;n) = T„ci*WgSP(-;n). 
// k is odd, then, for any sequence i = (ii, ^2, • • • , i2k), 

Proof. It is proved in a usual way. The Weingarten calculus for a Haar symplectic matrix 
5" gives 

n^h,^.<M■■■^^2,-u^J= E ^-(*) E Wg^^lr" V; n)f (t) 

where 

k 

P={Pl,-,P2k) '■=1 

for any r e S2k- □ 
The function Wg'~^^(-;?^) on S2k satisfies 

WgCi(C<;n) = e(C')WgCi(a;n) (a G ^2^1., CC e i^^), 

and hence Wg*-^^((;n) = for ( e Hj.. Furthermore, Wg^^(cr;n) = (a e S2k) if A; is 
odd. 

Example 5.2. 

Wg^'((}iii);n)=WgC^((}li^);n) 



2n+ 1 
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6 Conclusion 

We have made methods for computations of moments of matrix elements from classical 
compact Lie groups and classical compact symmetric spaces. Write C = A, B/D, C for 
unitary, orthogonal, symplectic groups, respectively. The moment for a classical group is 
given by the double sum 

(A-function in a) x (A-function in r) x Wg'' (cr^^'^r; n), 

cr T 

whereas that for a classical compact symmetric space is given by the single sum 

(A-function in cr) x Wg^{a; n). 

cr 

Here the A-function is 

• 5,(-, ■) defined in (E^D if C = A, A I, A II, A III, 

• A,(-) defined in (EJD if C = B/D, BD I, D III, 

• A;(-) defined in fl2:T2|) if C = C, C I, C II, 
and the Weingarten function Wg'' belongs to 

• Z{L{Sk)) if C = A, A III, 

. L{S2k,Hk) if C = B/D, A I, BD I, 

• L%S2k,Hk) ifC = C, All, CII. 

The Weingarten functions for C = D III or C I differ from others. In fact, if we let F{a) 
to be Wg°^^^((j; n) or Wg*"^((j"^; n), then F is the function on 5*2^ satisfying the property 

F{CO = e{OF{a) (a G ^2^, C, C e ^fc)- 

In particular, F identically vanishes on Hj^, and, moreover, on if k is odd. Weingarten 
functions except D III and C I have Fourier expansions in ^'^i oi' However, we 
could not find such expansions for these two cases. 
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